Comprehensive Approach to the Interpretation of the Electrical Properties of Film-Forming Molecules

This paper presents a general protocol for the interpretation of the electric surface potential of Langmuir monolayers based on a three-layer capacitor model. The measured values were correlated with the results from DFT molecular dynamics simulations, and, as a result, the local dielectric permittivities and dipole-moment components of molecules organized in the monolayer were obtained. The main advantage of the developed approach is applicability to amphiphiles of any type; irrespective of the structure of the polar head as well as the molecular organization and inclination in the surface film. The developed methodology was successively applied to an atypical surface-active compound, perfluorodecyldecane, and its derivatives containing the hydroxyl, thiol, and carboxyl moiety. The following contributions to the apparent dipole moments connected with the reorientation of water molecules and local dielectric permittivities in the vicinity of polar and apolar molecule parts, respectively, were determined: μw/εw = −0.85 D, εp = 5.00, and εa = 1.80. Moreover, the investigated perfluorodecyldecane derivatives were comprehensively characterized in terms of their surface activity, film rheology, and effective surface dissociation equilibria. The proposed methodology may be crucial for the process of the design and the preliminary characterization of molecules for sensor and material science applications.


INTRODUCTION
The classical method widely used to characterize monomolecular layers formed at the air/water interface (Langmuir monolayers) 1 is based on the measurement of surface pressure−area (π−A) isotherms. Analysis of the obtained π− A experimental curves provides first-hand information on the ability of the studied molecules to form stable monolayers, allows to determine (i) their characteristic parameters (such as molecular area), (ii) the monolayer physical state (based on the compressibility modulus C s −1 ), and (iii) collapse behavior. In addition, in the case of multicomponent monolayers, information on the mutual miscibility of film components and their interactions can be obtained. Apart from the classical surface manometry, many complementary techniques for in situ studies of monolayers were developed, involving microscopic 2 (fluorescence or polarizing microscopes), spectroscopic 3 (e.g., PM-IRRAS, SFG), and X-ray scattering 4 (e.g., GIXD) methods. The electric surface potential change (ΔV) measurements are equally important, although less frequently used. Although ΔV quantity meets the additivity rule and can be used to characterize miscibility and interactions between molecules in multicomponent systems, 5,6 it is mainly applied to investigate one-component monolayers. In this context, it enables the determination of the electrical parameters (such as dipole moments and electrical permittivity) of the molecules in the surface layer, 7 tracks changes in the molecular orientation during compression, and follows the formation of multilayers and domains. 8,9 Such information obtained for thin, organized, and defectless films is especially desired when considering applications of new compounds in biomaterials, nanoelectronics, and sensorics. 1,10 Although surface potential measurements have been performed for decades, theoretical models for relating experimental ΔV values to the ordering and structure of molecules in a film remain unchanged since the 1980s. Generally, an amphiphilic compound can be perceived as an electric dipole. Amphiphiles organized in the uniform monolayer at the air/water surface adopt characteristic orientation (the polar head is anchored in the water while the hydrophobic chains protrude to the air), which results in a charge gradient that is perpendicular to the surface. Thus, for the nonionized monolayer, the measured electric surface potential change ΔV is related to the normal (to the surface) component of the dipole moment of film molecules μ ⊥ where ε is the relative permittivity (dielectric constant) of the monolayer, ε 0 is the permittivity of vacuum, and A is the average area per molecule at the surface. The unknown ε value is often included in the so-called apparent dipole moment 11 of molecule in the film Theoretical models of surface potential (two-12 or threelayer 13,14 capacitors, reviewed in detail in refs 7, 15) are based on the Helmholtz equation, 16 where a monolayer is treated as a parallel-plate condenser comprising an array of uniformly distributed dipoles. In the most frequently used approach (based on the three-layer capacitor model; Demchak and Fort model 14 ), μ A can be expressed as a sum of three contributions In the above equation, ε p and ε a are local relative permittivities for polar and apolar parts of film molecules, respectively; μ ⊥ p and μ ⊥ a are the components of dipole moment normal to the water surface for polar and apolar groups, respectively; while μ ⊥ w /ε w is the contribution from the reorientation of water molecules. In the standard methodology, 17,18 surface electric potential change measurements are performed for the selected compounds bearing the same apolar parts and different polar parts, or vice versa. Then, the group dipole moments (μ ⊥ p , μ ⊥ a ) are calculated based on the structure of each amphiphile by adding up the bond dipole moments and taking into account the angles between them. The relative permittivity values (ε p , ε a ) are obtained by solving in pairs equations of type (3) for molecules of the same apolar parts and different polar parts, or vice versa. It is assumed that the contribution from the reorientation of water molecules for each pair of equations is the same. The described approach was improved by introducing DFT-optimized molecular conformations to determine μ ⊥ p and μ ⊥ a values; additionally, multiple linear regression was implemented to obtain ε p and ε a constants. 19 However, this methodology still has some limitations. First, it assumes the vertical orientation of the molecules at the air/ water interface in their most packed arrangement, while the molecules remain slightly inclined at this state in some cases. 20 Second, for nonclassical film-forming molecules with a nonamphiphilic structure, devoid of the typical polar group (for example, semifluorinated hydrocarbons), the differentiation between the polar and apolar parts of the molecule is problematic. Therefore, the aim of this research is to present a universal model that can be applied to molecules of any structure, based on DFT modeling. For our studies, we have selected a series of semifluorinated molecules, which have attracted a lot of attention in recent years due to their interesting properties and diverse applications. 21 These hybrid molecules comprise two incompatible fragments in their structure, hydrogenated and perfluorinated moieties, 20,22−26 and exhibit peculiar behavior. The most striking difference, compared to alkanes and perfluoroalkanes, is their surface activity 25−27 and liquid−crystalline properties, 25,26,28 which are absent in the counterpart apolar molecules. A plethora of papers on the physicochemical and structural properties of semifluorinated alkanes, including surface micelles formation, have appeared in the literature (for reviews, see 25,26,29 and references therein). However, the electrical properties of films of amphiphilic molecules containing a perfluorinated core attached to different polar groups have not been systematically studied. For our investigations, we have chosen 1,1,1,2,2,3, 3,4,4,5,5,6,6,7,7,8,8,9,9,10,10-henicosafluoroeicosane (F(CF 2 ) 10 (CH 2 ) 10 H abbr. F 10 H 10 ) and its amphiphilic derivatives containing different polar groups, namely, F 10 H 10 SH, F 10 H 10 OH, and F 10 H 10 COOH. The structures of the studied molecules are shown in Figure 1.

Langmuir Monolayer Characterization.
Surface pressure−area (π−A) and electrical surface potential change− area (ΔV−A) curves were measured simultaneously using NIMA equipment: a two-barrier trough of the total area of 600 cm 2 (612D) coupled with the surface pressure and surface potential sensor. In a typical experiment, 50−100 μL of the investigated compound solution in chloroform was carefully spread with a microsyringe onto the subphase surface. After solvent evaporation (approximately 10 min), the film was compressed with a barrier speed of 20 cm 2 /min. During the experiments, the subphase temperature was maintained at 20 ± 0.1°C by the Julabo thermostat. Surface pressure was recorded with an accuracy of ±0.1 mN/m using a Wilhelmy plate made of chromatography paper (Whatman Chr1) immersed in the subphase as the pressure probe. During electrical surface potential change measurements, the vibrating plate was located around 2 mm above the subphase surface while the reference electrode was placed in the subphase. The electrical surface potential was registered with an accuracy of ±15 mV and ±2 Å 2 /molecule. The surface pressure−area and electric surface potential change−area isotherms presented here are representative curves selected from at least two overlapping experiments.

Theoretical Calculations.
The dipole moments were calculated for previously geometrically optimized systems using the Gaussian 16 software package. 31 Geometry optimization was performed by density functional theory (DFT) modeling. All calculations were performed with the B3LYP functional, 32,33 a basis set including diffuse and polarization functions, i.e., 6-311++G(3df,3pd) 34,35 without damping. Systems were optimized with the default UltraFine integration grid, default integral cutoffs, and a combination of EDIIS and CDIIS tight convergence procedures, with no Fermi broadening. The dipole moments of polar μ ⊥ p and apolar μ ⊥ a parts of a molecule were determined utilizing the quantum theory of atoms in molecules (QTAIM) in AIMAll software. 36 Molecular dynamics calculations were performed in the Amber20 package. 37 Each analyzed system consisted of two symmetric rectangular monolayers, each having 128 perfluorodecyldecane or its derivative molecules, separated by 30,000 water molecules. Simulated systems were prepared in Packmol software. 38 The General AMBER Force Field 2 (GAFF 2) was used, with the partial atomic charge calculated by Gaussian 16 using the Hartree−Fock method and the 6-31G(d) basis set. Periodic boundary conditions were utilized. The TIP3P model 39 was used to simulate water molecules. The energy of the systems was minimized by 50,000 steps. The systems were equilibrated by 75,000 steps with a 0.001 ps timestep, followed by 300,000 steps with a 0.002 ps timestep. Production calculations were carried out in the isothermal−isobaric ensemble with the constant surface tension of 30 mN/m (NPγT) and with a 0.002 ps timestep. The temperature was set at 293 K and the Langevin thermostat was used. A Berendsen aerostat was used to control pressure at 1 bar. The simulation was carried out for 500 ns, and the last 10 ns were used for the analysis. Radial pair distribution functions were determined in the Cpptraj program. 40

Experimental π−A Isotherms.
In the first stage of our investigations, we looked at how the surface activity of the semifluorinated amphiphiles changes with the type of polar  Figure S1 (Supporting Information). As can be seen, the curves for F 10 H 10 , F 10 H 10 OH, and F 10 H 10 SH on water practically coincide with those recorded for the acidic subphase. In the case of F 10 H 10 COOH, the difference between isotherms is visible only at low surface pressure values (the lift off area of the isotherm recorded for the acidic subphase is ca. 2 Å greater compared to the isotherm registered for pure water). Therefore, it can be concluded that the surface activity of perfluorodecyldecane and its amphiphilic derivatives is hardly influenced by the degree of dissociation. To analyze how the surface activity of perfluorodecyldecane derivatives varies with the type of polar group, the representative isotherms recorded on the HCl solution together with the calculated compressibility moduli (C s −1 ) curves have been compared in Figure 2. The characteristic parameters read from the plots are summarized in Table 1. Values of C s −1 , calculated on the basis of the experimental isotherm datapoints (applying , are helpful in describing the physical state of the surface film. Namely, C s −1 values below 25 mN/m suggest that the film is in a low-density liquid phase; the ranges of 25−50 and 100−250 mN/m are characteristics of the liquid expanded and liquid condensed states, respectively, while the film is in the solid state for C s −1 above 500 mN/m. 1 The π−A isotherm of F 10 H 10 starts to rise at the molecular area below 27.8 Å 2 and collapses at an area of 19.0 Å 2 / molecule and a surface pressure of 9.5 mN/m. The isotherm shows a clear inflection above the surface pressure of 3 mN/m, which is also clearly visible in the C s −1 −A plot (as a minimum), and the calculated compressibility moduli values suggest that the inflection is due to the phase transition within liquid state. 11,24 The introduction of a hydrophilic group in the terminal position of the perfluorodecyldecane moiety strikingly influences the surface activity of the compounds obtained. Namely, compared to F 10 H 10 , isotherms of its amphiphilic derivatives show some common features: (i) the lift off areas are shifted to larger values (above 41 Å 2 /molecule), (ii) the slope of the curves is more vertical and without visible inflections (the monolayers are in a liquid condensed state), (iii) the collapse pressure is significantly higher. Nevertheless, the analysis of the π−A curves also shows some differences between the curves recorded for functionalized perfluorodecyldecanes. F 10 H 10 SH forms slightly more expanded films compared to F 10 H 10 OH and F 10 H 10 COOH. This is because the measured π−A dependencies for F 10 H 10 SH start to increase at larger areas (approximately 41.9 Å 2 /molecule) than for F 10 H 10 OH and the slope of the isotherm is more inclined than for F 10 H 10 OH and F 10 H 10 COOH. It is also reflected in the calculated compressional moduli values, which are almost 1.7 times smaller than those for F 10 H 10 OH and F 10 H 10 COOH; however, the physical state of the monolayers remains the same (liquid condensed). Furthermore, the collapse pressure of F 10 H 10 SH film is quite low (approximately 25.4 mN/m). The F 10 H 10 OH and F 10 H 10 COOH isotherms have a similar slope (and similar maximum values of C s −1 ); however, the collapse pressure value of F 10 H 10 OH is higher than that of F 10 H 10 COOH (59.8 vs 50.0 mN/m).

Experimental ΔV−A Isotherms.
In the next step, we examined the electrical properties of surface films formed by perfluorodecyldecane and its derivatives at the air−water boundary. Initially, the electric surface potential change−area (ΔV−A) isotherms were measured for films formed on the ultrapure water as well as on 0.001 mol/dm 3 aqueous HCl solution ( Figure S2, Supporting Information). It was noticed that the experimental curves registered on both subphases overlap only for F 10 H 10 . However, in the case of other compounds (F 10 H 10 SH, F 10 H 10 OH, and F 10 H 10 COOH), there is an additional contribution to the electric surface potential connected with the double-layer potential (ψ 0 ). According to the Gouy−Chapman theory, the double-layer potential (ψ 0 ) is a result of the dissociation of polar head groups at the water− air interphase and/or their involvement in a hydrogen-bonding network. 42 This issue is discussed in detail in Section 3.4.
To determine the electric properties of film molecules resulting directly from their structure, the electric surface potential change−area (ΔV−A) curves measured for monolayers formed on 0.001 mol/dm 3 aqueous HCl solution were used (Figure 3). In this way, the double-layer contribution (ψ 0 ) resulting from different surface ionization degrees of polar groups in the investigated compounds can be avoided. The characteristic parameters read from the experimental curves are summarized in Table 1.
The ΔV−A dependencies measured for all investigated compounds are characterized by a similar course. At the beginning of compression, the ΔV values remain constant at approximately zero until the so-called critical area (A c ) 7,43 is reached. The value of A c is different for each of the compounds tested and indicates the point when the hydrogen bonds with the water molecules are broken and the monolayer begins to The value of the area per molecule corresponding to the π−A isotherm's lift off (lift off point). b The value of the area per molecule corresponding to the beginning of the decrease in the ΔV−A isotherm (critical area). c The maximum values of the compressibility moduli. d The area per molecule corresponding to the maximum value of the compressibility moduli. e The electric surface potential corresponding to the maximum packing of monolayer molecules.
The Journal of Physical Chemistry B pubs.acs.org/JPCB Article organize. 7 As expected, A c has the smallest value for a purely hydrophobic molecule, F 10 H 10 , which is not involved in hydrogen bonding with water molecules. After the value of A c upon compression is exceeded, the electric surface potential isotherm gradually decreases toward more negative values until it reaches a clear inflection point at small areas per molecule. It is interesting that for F 10 H 10 , the surface potential isotherm raises toward less negative values upon further compression. This observation, in addition to the inflection in the π−A isotherm, proves that the F 10 H 10 molecules in the monolayer undergo a phase transition during compression. 11 In addition to the critical area (A c ), another key parameter that can be obtained from the experimental ΔV−A curves is ΔV max value, which is the electric surface potential change corresponding to the closest packing in the monolayer. ΔV max value is read from the ΔV−A curve for the molecular area A max , which corresponds to the maximum of compressibility moduli. Values of experimental apparent dipole moments (μ A exp ) can be calculated from experimental ΔV−A dependencies by applying the following equation where A max and ΔV max are the area per molecule and the electric surface potential corresponding to the maximum packing of the monolayer molecules. These values obtained for the studied semifluorinated molecules are summarized in Table 1.
As can be seen, the apparent dipole moment for F 10 H 10 is small and differs from the values for amphiphilic derivatives (F 10 H 10 SH, F 10 H 10 OH, F 10 H 10 COOH). This can be explained by taking into account the molecular structure of the investigated compounds, which is discussed in the next section.

Interpretation of Experimental Apparent Dipole Moment Values of Monolayers Formed by Nonionized
Amphiphiles. Taking into account the contributions to the apparent dipole moment of the film molecule (eq 3), the normal components of the dipole moment of the nonionized polar and apolar parts of the molecule (μ ⊥ p and μ ⊥ a ) can be calculated based on its structure. The standard approach involves optimization of the molecular conformation in a vacuum using semiempirical or DFT methods. The values of μ ⊥ p and μ ⊥ a for the polar and apolar parts of the molecule can be determined assuming that the investigated amphiphile adopts vertical orientation in a closely packed monolayer at the water−air interphase. Such an assumption was successfully applied to aromatic carboxylic acids 44 and phosphocholines. 19 However, for perfluorodecyldecane derivatives, this approach failed possibly because of two main reasons. First, the investigated molecules were found to be inclined in a closely packed film at the surface of the water, as shown in ref 11. Therefore, the angle between the main axis of the molecule (defined as the vector that connects the extreme positions of carbon atoms in an aliphatic chain) in the closely packed arrangement and the normal to the water−air interphase (θ) should be taken into account. Second, the precise distinction between the polar and apolar parts of the molecule may be problematic; therefore, information on the hydration shell of the hydrophilic groups is desirable. To investigate these issues, the simulations of molecular dynamics were carried out for the investigated perfluorodecyldecane derivatives with the assumption of nonionized polar groups. The obtained results confirmed that the orientation of the investigated compounds in their most packed arrangements (corresponding to max C s −1 ) is not vertical (Figure 4); the deviation of the main axis of the molecules from the normal to the water−air interphase (θ) is approximately 53°for F 10 H 10 SH, while that for the remaining studied molecules is below 47°( Table 2).
To identify the hydrated parts of the molecules simulated in molecular dynamics, we determined the radial pair distribution functions between carbon and oxygen atoms in the investigated molecules and oxygen atoms in water molecules ( Figure 5). It was assumed that the maximum of the radial distribution function corresponding to a distance slightly above   The Journal of Physical Chemistry B pubs.acs.org/JPCB Article 3 Å suggests the presence of the selected atom in the first adsorbed water layer. Otherwise, a value greater than 4 Å suggests that the neighboring atom and/or further atoms in the molecule are hydrated. As it can be noticed, for each investigated compound, the distance between the first carbon atom in the hydrocarbon chain of the amphiphile and the oxygen atoms of the water molecules is clearly defined and equals approximately 3 Å. This suggests that the following groups can be perceived as hydrated (polar): −CH 3 (in F 10 H 10 ), −CH 2 SH (in F 10 H 10 SH), −CH 2 OH (in F 10 H 10 OH), and −COOH (in F 10 H 10 COOH). Taking into account the inclination of the investigated amphiphiles at the surface in their most packed arrangements and the assignments of the hydrated parts of the molecules, the values of μ ⊥ p and μ ⊥ a were determined using the DFT approach (see Table 2).  , ε p = 5.00, and ε a = 1.80. The calculated local relative permittivities are slightly lower in comparison to the literature values, which lie in the range 6−7 and 2−3 for ε p and ε a , respectively. 45 Nonetheless, the calculated local dielectric permittivity values prove a good correlation between the experimental values of apparent dipole moments and those obtained from DFT molecular dynamics.

Interpretation of Experimental Apparent Dipole Moments of Monolayers from Ionized Amphiphiles.
In the next stage of our investigations, the electrical properties resulting exclusively from polar interactions and/or H-bonding of the polar group attached to the perfluorodecyldecane moiety were analyzed. For this purpose, the ΔV−A isotherms on subphases differing in ionic strength (containing the following concentrations of NaCl: 0.1, 0.01, and 0.001 mol/dm 3 ) were measured for each investigated compound. The resulting ΔV− A curves were compared with isotherms registered using an aqueous subphase containing HCl (0.001 mol/dm 3 ) and are presented in Figure 6. Additionally, the ΔV max values (corresponding to the closest monolayer packing) were read and are compiled in Table 3.
As it can be seen, ΔV max values for all investigated compounds follow the same trend depending on the subphase applied. The ΔV max values for the subphase containing 0.001 mol/dm 3 HCl are the highest (the least negative) and correspond to the uncharged monolayer. The surface potential change for charged monolayers is the lowest (most negative) for the subphase containing 0.001 mol/dm 3 NaCl and increases gradually as the subphase ionic strength increases. As already mentioned, for the partially charged surface layer, the electric surface potential change (ΔV i ) is not only related to the molecular dipole moments but also the double-layer potential (ψ 0 ) should be taken into consideration Thus, knowing that the first segment in eq 6 corresponds to the surface potential of nonionized monolayer, ψ 0 can be calculated where ΔV max i and ΔV max n are the electric surface potential change values corresponding to the closest packing in the ionized and nonionized monolayer, respectively. The obtained ψ 0 values are summarized in Table 3, and their dependence on the ionic strength of the subphase is plotted in Figure 7.
From the initial analysis, it should be noticed that ψ 0 values calculated for F 10 H 10 SH are the lowest and only one of them surpasses the experimental inaccuracy value. Nonetheless, the analysis of the data in Figure 7 shows that ψ 0 values for all functionalized perfuorodecyldecane derivatives decrease with the increasing ionic strength. For each of the investigated compounds, the slope of this trend is different. To better understand these differences, the ionization degrees and dissociation constants of molecules at the surface were calculated using the literature methodology. 42,46 The surface ionization degree at a temperature of 20°C can be obtained from the following equation where f NaCl is activity coefficient value equal to 1.00, 0.90, or 0.79 for NaCl concentration of 0.001, 0.01, or 0.1 mol/dm 3 , respectively. Finally, the pK a S values were calculated using the Henderson−Hasselbalch equation   The knowledge of pK a S value is of utmost importance especially in biomedical sciences and sensorics because the majority of pivotal processes take place at the interfaces, where the value of the dissociation constant for an ionizable molecule differs from that in the bulk phase. For aliphatic (unsubstituted) carboxylic acids, bulk pK a falls within the range of 4.7−5.0, 47 while pK a S at the air/water interface is higher. For example, the value pK a S = 5.7 ± 0.2 was determined for condensed stearic acid in Langmuir monolayers from surface potential measurements. 42,48 Similarly, for nonanoic acid, it is equal to 5.9 ± 0.6 as obtained from a combination of surface tension measurements and molecular dynamics simulations. 49 However, it is worth to mention that in the case of amphiphiles that are slightly soluble in water (i.e., nonanoic acid), the surface adsorption/desorption equilibrium should also be taken into consideration. 49,50 For aliphatic bases in bulk, the value reads 10.6−10.9, 47 while it is 10.1 at the interface (determined for nonadecylamine). 51 Similar results were obtained applying different methods, e.g., SFG. 52,53 It is clear that at the interface there is a shift of pK of about 0.7 for both acids and bases toward neutral pH. This is due to the fact that the conditions at the phase boundary are more favorable for a larger accumulation of uncharged molecules, while ionized molecules avoid contact with the interface. It can be noticed that pK a S values determined for amphiphilic derivatives of perfluorodecyldecane are noticeably lower than for analogous hydrocarbon derivatives, i.e., the average pK a S for F 10 H 10 COOH is ca. 5.2, while that for stearic acid is 5.7. This can be explained by the electron-withdrawing effect of fluorine atoms, which causes the pK a for fluorinated compounds to be significantly lower as compared to those of their aliphatic analogues. For example, the substitution of one hydrogen atom in acetic acid by a fluorine atom decreases pK a from 4.8 to 2.6 in bulk. 47 Similarly, the pK a value for ethanol is equal to ca. 16.0, while for 2,2,2trifluoroethanol, it drops to 12.4. 54 In our case, the electronwithdrawing effect of fluorine atoms is less pronounced due to the quite long distance between the perfluorinated segment and the ionizable group in the molecule. Furthermore, the comparison of the determined pK a S values demonstrates that the ability of the investigated molecules to undergo dissociation at the water/air boundary increases in the order F 10 H 10 ≪ F 10 H 10 SH < F 10 H 10 OH ≪ F 10 H 10 COOH. Our results show that effective surface ionization constants for the investigated fluorinated thiol and alcohol are comparable, whereas in bulk, thiols have been reported to be stronger electrolytes than alcohols (for example, bulk pK a for ethanol and ethanethiol are equal to 16.0 54 and 10.9, 55 respectively); however, this relation has been obtained for hydrogenated molecules.

CONCLUSIONS
Electric surface potential measurements are much less frequently applied for the classical characterization of amphiphiles in monolayers at the free water surface, which is mainly due to difficulties in the interpretation of the measured ΔV values. In this paper, a general protocol for estimating particular contributions to a three-layer capacitor model 14 was developed based on a correlation of the experimentally obtained values with the results from DFT molecular dynamics simulations. The presented methodology has some undeniable advantages in comparison to the previous approaches (for a detailed discussion of previous approaches, see Section 1). First, it can be applied to every amphiphile (even to the surface-active compounds that are devoid of a typical polar group(s) in their structure) as radial distribution functions provide information on the hydration shell. For example, for typical amphiphiles, we found that the hydrated polar part can be considered as −CH 2 OH and not −OH, or −CH 2 SH, and not −SH. Simultaneously, the methyl group in perfluorodecyldecane can also be considered to be hydrated. Second, no assumption of the vertical orientation of molecules at the surface is necessary as molecular dynamics simulations give an average angle between the molecule axis and normal to the water surface. As a result, it becomes possible to correlate the experimental results with the organization, inclination, and average conformation of molecules of interest in the monolayer. Moreover, for ionizable compounds, it has been demonstrated that information on surface dissociation equilibrium is also available. The above-mentioned benefits of combining electrical surface potential measurements with molecular dynamics simulations may be of key importance for the design and initial characterization of molecules for sensor and materials science applications.
Additional figures presenting the surface pressure−area and electric surface potential change−area isotherms (PDF)